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c3 ■ By combining high-resolution HST and wide-field ground based observations, 

in ultraviolet and optical bands, we study the Blue Stragglers Star (BSS) popula- 
tion of the galactic globular cluster M5 (NGC 5904) from its very central regions 
up to its periphery The BSS distribution is highly peaked in the cluster center, 
decreases at intermediate radii and rises again outward. Such a bimodal dis- 
tribution is similar to those previously observed in other globular clusters (M3, 
47 Tucanae, NGC 6752). As for these clusters, dynamical simulations suggest 
that, while the majority of BSS in M5 could be originated by stellar collisions, 
a significant fraction (20-40%) of BSS generated by mass transfer processes in 
primordial binaries is required to reproduce the observed radial distribution. A 
candidate BSS has been detected beyond the cluster tidal radius. If confirmed, 
this could represent an interesting case of an "evaporating" BSS. 
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Subject headings: Globular clusters: individual (M5); stars: evolution - binaries: 
general - blue stragglers 



1. INTRODUCTION 

In globular cluster (GC) color-magnitude diagrams (CMD) blue straggler stars (BSS) 
appear to be brighter and bluer than the Turn-Off (TO) stars and lie along an extension of 
the Main Sequence. Since BSS mimic a rejuvenated stellar population with masses larger 
than the normal cluster stars (this is also confirmed by direct mass measurements; e.g. Shara 
et al. 1997), they are thought to be objects that have increased their initial mass during 
their evolution by means of some process. Two main scenarios have been proposed for their 
formation: the collisional scenario suggests that BSS are the end-products of stellar mergers 
induced by collisions (COL-BSS), while in the mass-transfer scenario BSS form by the mass- 
transfer activity between two companions in a binary system (MT-BSS), possibly up to the 
complete coalescence of the two stars. Hence, understanding the origin of BSS in stellar 
clusters provides valuable insight both on the binary evolution processes and on the effects 
of dynamical interactions on the (otherwise normal) stellar evolution. 

The re lative efficiency of the two formation mechan i sms is thought to depend on the en- 



viron ment (IFusi Pecci et al. Ill992l ; iFerraro et al.lll999al ; iBellazzini et al.ll2002l ; iFerraro et al. 



20031 ). COL-BSS are expected to be formed preferentially in high-density environments (i.e., 
the GC central regions), where stellar collisions are most probable, and MT-BSS should 
mainly populate lower density environments (the cluster peripheries), where binary systems 
can more easily evolve in isolation without suffering exchanges or ionization due to gravita- 
tional encounters. The overall scenario is complicated by the fact that primordial binaries 
can also sink to the core due to mass segregation processes, and "new" binaries can be formed 
in the cluster centers by gravitational encounters. The two formation mechanisms are likely 
to be at work simultaneously in every GC (see the case of M3 as an example; Ferraro et al. 
1993, 1997), but the identification of the cluster properties that mainly affect their relative 
efficiency is still an open issue. 

One possibility for distinguishing between the two types of BSS is offered by high- 
resolution spectroscopic studies. Anomalous chemical abundances are expected at the surface 



1 Based on observations with the NASA/ESA HST, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS5-26555. Also based on WFI observations 
collected at the European Southern Observatory, La Silla, Chile, within the observing programs 62.L-0354 
and 64.L-0439. 
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of BSS resulting from MT activi ty (Sarna fc de Greve 1996h , whil e they are not predicted 
in case of a collisional formation (jLombardi. Rasio &: Shapiro! 1995 ). Such stu dies have just 



becom e feasible, and the results found in the case of 47 Tucanae (47 Tuc; iFerraro et al. 



2006al ) are encouraging. The detection of unexpected properties of stars along standard 
evolutionary sequences (e.g., variability, anomalous population fractions, or peculi ar radial 



distributions) can help estimating the fraction of binaries within a cluster (see, e.g.. iBailyn 



1994 ; lAlbrow et al.ll200ll ; iBellazzini et al.l 120021 ; iBeccari et al.ll2006l ). but such evidence does 
not directly allow the determination of the relative efficiency of the two BSS formation 
processes. 

The most widely a pplicable tool to probe the origin of BSS is their radial distribution 



within the clusters (see IFerraro 1120061 . for a review). This has been observed to be bimodal 



(i.e., highly peaked in the cluster cente rs and peripheries, a nd signific antly lower at inter - 
mediate ra dii) in at least 4 GCs: M3 ( Ferraro et al.lll997l). 47 Tuc (jFerraro et al.l 12004 ). 
NGC 6752 JSabbi et al.ll2004j ). and M5 JWarren. Sandauist fe Bolte Ibood. hereafter W06) 
Preli minary evidence of bimodality has also been found in M55 (jZaggia. Piotto fc Capaccioli 
19971 ). Dynamical simul ations suggest that t he bi modal radial distributions observed in M3, 
47 Tuc and NGC 6752 JMapelli et al.lbooi l2006h result from ~ 40 - 50% of MT-BSS with 
the balance being COL-BSS. In this context, the case of ui Cen is atypical: the BSS radial 
distribution in this cluster is flat (IFerraro et al.ll2006bl ). and mass segregation processes have 
not yet play ed a major role, thu s implying that this system is populated by a vast majority 
of MT-BSS (jMapelli et al.ll2006l ). These results demonstrate that detailed studies of the BSS 
radial distribution within GCs are very powerful tools for better understanding the complex 
interplay between dynamics and stellar evolution in dense stellar systems. 

In the present paper we extend this kind of investigation to M5 (NGC 5904). With HST- 
WFPC2 and -ACS ultraviolet and optical high-resolution images of the core we have been 
able to efficiently detect the BSS population even in the severely crowded central regions. 
Moreover, with wide-field optical observations performed with ESO-WFI we sampled the 
entire cluster extension. The combination of these two data sets allowed us to study the 
dynamical properties of M5, accurately redetermining its center of gravity, its surface density 
profile, and the BSS radial distribution over the entire cluster. The BSS population of M5 
has been recently studied by W06, but we have extended the analysis to larger distances 
from the cluster center, and we have used Monte-Carlo dynamical simulations to interpret 
the observational results. 
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2. OBSERVATIONS AND DATA ANALYSIS 
2.1. The data sets 

The present study is based on a combination of two different photometric data sets: 

1. The high-resolution set - It consists of a series of ultraviolet (UV) and optical images 
of the cluster center obtained with HST-WFPC2 (Prop. 6607, P.I. Ferraro). To efficiently 
resolve the stars in the highly crowded central regions, the Planetary Camera (PC, being the 
highest resolution instrument: O'.'046/pixel) has been pointed approximately on the cluster 
center, while the three Wide Field Cameras (WF, having a lower resolution: O'.'l/pixel) have 
been used to sample the surrounding regions. Observations have been performed through 
filter F255W (medium UV) in order to efficiently select the BSS and horizontal branch 
(HB) populations, and through filters F336W (approximately corresponding to an U filter) 
and F555W (V) for the red giant branch (RGB) population and to guarantee a proper 
combination with the ground-based data set (see below). The photometric reduction of the 
high-resolution images was carried out using ROMAFOT (Buonanno et al. 1983), a package 
developed to perform accurate photometry in crowded fields and specifically optimized to 
handle under-sampled Point Spread Functions (PSFs; Buonanno & Iannicola 1989), as in 
the case of the HST-WF chips. 

To obtain a better coverage of the innermost regions of the cluster, we have also used a 
set of public HST-WFPC2 and HST-ACS observations. T he HST-WFPC2 dat a set has been 



obtained through filters F439W (B) and F555W (V) bv IPiotto et al. I (120021 ) . and because 
of the different orientation of the camera, it is complementary to ours. Additional HST-ACS 
data in filters F435W (B), F606W (V), and F814W (/) have been retrieved from the ESO- 
STECF Science Archive, and have been used to sample the central area not covered by the 
WFPC2 observations. All the ACS images were properly corrected fo r geometric dis t ortion s 



and effective flux (over the pixel area) following the prescriptions of ISirianni et al.l (120051 ). 



The photometric analysis was performed independently in the three drizzled images by using 



the aperture photometry code SExtractor (Source- Extractor, iBertin fc Arnoutd 119961 ). and 
adopting a fixed aperture radius of 2.5 pixels (0.125"). The magnitude lists were finally 
cross-correlated in order to obtain a combined catalog. The adopted combination of the 
three HST data sets is sketched in Figure CD and provided a good coverage of the cluster up 
to r = 115". 

2. The wide-field set - A complementary set of wide-field B and V images was secured by 
using the Wide Field Imager (WFI) at the 2.2m ESO-MPI telescope during an observing run 
in April 2000. Thanks to the exceptional imaging capabilities of WFI (each image consists of 
a mosaic of 8 CCDs, for a global field of view of 34' x 34'), these data cover the entire cluster 
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extension (see Figured where the cluster is roughly centered on CCD #7). The raw WFI 
images were corrected for bias and flat field, and the overscan regions were trimmed using 
IRAF@ t ools. The PSF fitting procedure was performed independently on each image using 
DoPhot (jSchechter. Mateo fc Saha" 1993 ). All the uncertain detections, usually caused by 
photometric blends, stars near the CCD gaps or saturated stars, have been checked one by 
one using ROMAFOT (Buonanno et al. 1983). 



2.2. Astrometry and center of gravity 

The HST+WFI catalog has been placed on the absolute astrometric system by adopting 
the procedure already described in Ferraro et al. (2001, 2003). The new astrometric Guide 
Star Catalog (GSC-IH) was used to search for astrometric standard stars in the WFI field of 
view (FoV), and a cross-correlation tool specifically developed at the Bologna Observatory 
(Montegriffo et al. 2003, private communication) has been employed to obtain an astrometric 
solution for each of the 8 CCDs. Several hundred GSC-II reference stars were found in each 
chip, thus allowing an accurate absolute positioning of the stars. Then, a few hundred stars 
in common between the WFI and the HST FoVs have been used as secondary standards to 
place the HST catalog on the same absolute astrometric system. At the end of the procedure 
the global uncertainties in the astrometric solution are of the order of ~ 0'/2, both in right 
ascension (a) and declination (5). 

Given the absolute positions of individual stars in the innermost regions of the cluster, 
the center of gravity C grav has been determined by averaging coordinates a and 5 of all 
stars lying in the PC FoV following the iterative procedure described in Montegriffo et 
al. (1995; see also Ferraro et al. 2003, 2004). In order to correct for spurious effects 
due to incompleteness in the very inner regions of the cluster, we considered two samples 
with different limiting magnitudes (m 555 < 19.5 and m 555 < 20), and we computed the 
barycenter of stars for each sample. The two estimates agree within ~ 1", giving C grav at 
a(J2000) = 15 h 18 m 33!53, 5(J2000) = +2° 4' 57^06, with a lex uncertainty of 0"5 in both a 
and 5, corresponding to about 10 pixels in the PC image. This value of C grav is l ocated at 



~ 4" south-west (Aa = —4", A5 = — 0'.'9) from that previously derived by lHarris I (119961 ) on 
the basis of the surface brightness distribution. 



2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



3 Available at 



http : //www-gsss . stsci . edu/Catalogs/GSC/GSC2/GSC2 . htm 
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2.3. Photometric calibration and definition of the catalogs 



The optical HST magnitudes (i.e., those obtained through the WFPC2 filters F439W 



Dolphin I (120001 ) . 



and F555W, and through ACS filters F435W, F606W, F 814W), as well as 
V magnitudes have been all calibrated on the catalog of ISandauist et a 
magnitud es m 1fin and ha ve been calibrated to the iHoltzman et al 
following 



the WFI B and 



1996h . The UV 
"fll995h zero-points 



erraro et al.l (119971 . 2001), while the U magnitude m 336 has been calibrated to 



In order to reduce spurious effects due to the low resolution of the ground-based obser- 
vations in the most crowded regions of the cluster, we use only the HST data for the inner 
115", this value being imposed by the FoV of the WFPC2 and ACS cameras (see Figure 
[p. In particular, we define as HST sample the ensemble of all the stars in the WFPC2 
and ACS combined catalog having r < 115" from the center, and as WFI sample all stars 
detected with WFI at r > 115" (see Figure [2]). The CMDs of the HST and WFI samples in 
the (V, U — V) and (V, B — V) planes are shown in Figure [3J 



2.4. Density profile 

We have determined the projected density profile over the entire cluster extension, from 
Cgrav out to ~ 1400" ~ 23'3, by direct star counts, considering only stars brighter than 
V = 20 (see Figure [3]) in order to avoid incompleteness biases. The brightest RGB stars that 
are strongly saturated in the ACS data set have been excluded from the analysis, but since 
they are few in number, the effect on the resulting density profile is completely negligible. 
Following the procedure already described in Ferraro et al. (1999a, 2004), we have divided 
the entire HST+WFI sample in 27 concentric annuli, each centered on C grav and split in 
an adequate number of sub-sectors. The number of stars lying within each sub-sector was 
counted, and the star density was obtained by dividing these values by the corresponding 
sub-sector areas. The stellar density in each annulus was then obtained as the average of 
the sub-sector densities, and its standard deviation was estimated from the variance among 
the sub-sectors. 

The radial density profile thus derived is plotted in FigureHJ where we also show the best- 
fit mono-mass King model and the corresponding values of the core radius and concentration: 
r c = 27" (with a typical error of ~ ±2") and c = 1.68, respectively. These values confirm 
that M5 h as not yet experienced core co ll apse, and they are in good agreement with those 



consistent with those listed by 



Harris 



quoted by iMcLaughlin fc van der Marel (120051 . r c = 26" 3 and c = 1.71), and marginally 



19961 . r c = 25" 2 and c = 1.83), both derived from 
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the surface brightness profile. Our value o f ry corresponds to y 1 pc assuming the distance 



modulus (m — M) = 14.37 (d ~ 7.5 Kpc, iFerraro et al.lll999bl ). 



3. DEFINITION OF THE SAMPLES 

In order to study the BSS radial distribution and detect possible peculiarities, both the 
BSS and a reference population must be properly defined. Since the HST and the WFI data 
sets have been observed in different photometric bands, different selection boxes are needed 
to separate the samples in the CMDs. The adopted strategy is described in the following 
sections (see also Ferraro et al. 2004 for a detailed discussion of this issue). 



3.1. The BSS selection 

At UV wavelengths BSS are among the brightest objects in a GC, and RGB stars are 
particularly faint. By combining these advantages with the high-resolution capability of HST, 
the usual problems associated with photometric blends and crowding in the high density 
central regions of GCs are minimized, and BSS can be most reliably recognized and separated 
from the other populations in the UV CMDs. For these reasons our primary criterion for the 
definition of the BSS sample is based on the position of stars in the (m 255 , m 255 — U) plane. 
In order to avoid incompleteness bias and the possible contamination from TO and sub-giant 
branch stars, we have adopted a limiting magnitude m 255 = 18.35, roughly corresponding 
to 1 magnitude brighter than the cluster TO. This is also the limiting magnitude used by 
W06, facilitating the comparison with their study. The resulting BSS selection box in the 
UV CMD is shown in Figure O Once selected in the UV CMD, the bulk of the BSS lying in 
the field in common with the optical-HST sample has been used to define the selection box 
and the limiting magnitude in the (B, B — V) plane. The latter turns out to be B ~ 17.85, 
and the adopted BSS selection box in the optical CMD is shown in Figure [6j The two stars 
lying outside the selection box (namely BSS-19 and BSS-20 in Table [TJ have been identified 
as BSS from the (m 2 55, m 255 — U) CMD. Indeed, they are typical examples of how the optical 
magnitudes are prone to blend/crowding problems, while the BSS selection in UV bands is 
much more secure and reliable. An additional BSS (BSS-47 in Table [T]) lies near the edge 
of the ACS FoV and has only V and I observations; thus it was selected in the (V, V — I) 
plane (see Figured where this BSS is shown together with the other 5 identified in the ACS 
complementary sample). 

With these criteria we have identified 60 BSS: 47 BSS in the HST sample (r < 115") 
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and 13 in the WFI one. Their coordinates and magnitudes are listed in Table [TJ Out of the 
47 BSS identified in the HST sample, 41 are from the WFPC2 data set, and 6 from the ACS 
catalog. As shown in Figure [T] their projected distribution is quite asymmetric with the N-E 
sector seemingly underpopulated. The statistical significance of such an asymmetry appears 
even higher if only the BSS outside the core are considered. However a quantitative discussion 
of this topic is not warranted unless additional evidences supporting this anomalous spatial 
distribution are collected. One of the inner BSS (BSS-29 in Tabled]) lying at 21 / /76 from the 



center , corresponds to the low-amplitude variable HST-V28 identified by iDrissen fc Shara 
fll998H 4 l. In the WFI sample (r > 115") we find 13 BSS, with a more symmetric spatial 



distribution (see Figure [2]). The most distant BSS (BSS-60 in Table [fl marked with an 
empty triangle in Figj6]) lies at ~ 24' from the center, i.e., beyond the cluster tidal radius. 
Hence, it might be an evaporating BSS previously belonging to the cluster. However, further 
investigations are needed before firmly assessing this issue. 

In order to perform a proper comparison with W06 study, we have transformed their 
BSS catalog in our astrometric system, and we have found that 50 BSS of their bright sample 
lie at r < 115": 35 are from the HST sample, 13 from the Canada France Hawaii Telescope 
(CFHT) data set, and 2 from the Cerro Tololo Inter-American Observatory (CTIO) sample; 
in the outer regions (115" < r ^ 425") 9 BSS are identified, all from the CTIO data set. 

By cross correlating W06 bright sample with our catalog we have found 43 BSS in 
common (see Tabled]), 37 at r < 115" and 6 outward. In particular, 33 BSS out of the 41 
(i.e., 80% of the total) that we have identified in the WFPC2-HST sampl^f] are found in 
both catalogs, while 3 of our BSS belong 5 to their faint BSS sample (namely, BSS-27, 34, 
and 40, corresponding to their Core BSS 70, 79, and 76, respectively), 5 of our BSS have 
been missed in W06 paper, and 2 objects in their sample are classified as HB stars in our 
study. This is probably due to different selection criteria, and/or small differences in the 
measured magnitudes, caused by the different data reduction procedures and photometric 
analysis. For example, W06 identify the BSS on the basis of both the UV and the optical 
observations, while we select the BSS only in the UV plane whenever possible. Out of the 
other 15 BSS found at r < 115" in the ground-based CFHT/CTIO sample of W06, 8 BSS 
(Core BSS 38-45 in their Table 2) clearly are false identifications. They are arranged in a 
very unlikely ring around a strongly saturated star, as can be seen in Figure [8] where the 
position in the sky of the 8 spurious BSS are overplotted on the CFHT image. Though they 



The observations presented here do not have the time coverage needed to properly search for BSS 
variability. 

5 Note that the WFPC2-HST observations used in W06 and in the present study are the same. 
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clearly are spurious identifications, they still define a clean sequence in the (B, B — I) CMD, 
nicely mimicking the BSS magnitudes and colors. As already discussed in previous papers, 
this once again demonstrates how automatic procedures for the search of peculiar objects 
are prone to errors, especially when using ground-based observations to probe very crowded 
stellar regions. We emphasize that all the candidate BSS listed in our Table [I] have been 
visually inspected evaluating the quality and the precision of the PSF fitting. This procedure 
significantly reduces the possibility of introducing spurious objects in the sample. Out of the 
remaining 7 BSS, 4 objects (namely their Core BSS 32, 30, 37 and 28) are also confirmed 
by our ACS observations (BSS-42, 43, 44, and 45 respectively), while 2 others (their Core 
BSS 27 and Ground BSS 6) are not found in the ACS data set, and the remaining one (their 
Ground BSS 7) is not included in our observation FoV. In turn, two BSS identified in our 
ACS data set (BSS 46 and 47) are missed in their sample. Concerning the BSS lying at 
115" < r < 450", 6 objects (out of 9 found in both samples) are in common between the two 
catalogs (see Tabled]), one (BSS-55) belongs to W06 faint sample (their Ground BSS 23), 
while the remaining 2 do not coincide. Moreover, 4 additional BSS have been identified at 
r > 450" in our study. 



3.2. The reference population 

Since the HB sequence is bright and well separable in the UV and optical CMDs, we 
chose these stars as the primary representative population of normal cluster stars to be used 
for the comparison with the BSS data set. As with the BSS, the HB sample was first defined 
in the (771255, m 255 — U) plane, and the corresponding selection box in (B, B — V) has then 
been determined by using the stars in common between the UV and the optical samples. The 
resulting selection boxes in both diagrams are shown in Figures [5] and [6J and are designed 
to include the bulk of HB stara§. Slightly different selection boxes would include or exclude 
a few stars only without affecting the results. 

We have used WFI observations to roughly estimate the impact of possible foreground 
field stars contamination on the cluster population selection. As shown in the right-hand 
panel of Figure [6], field stars appear to define an almost vertical sequence at 0.4 < B — V < 1 
in the (B, B — V) CMD. Hence, they do not affect the BSS selection box, but marginally 
contaminate the reddest end of the HB. In particular, 5 objects have been found to lie within 
the adopted HB box in the region at r > r t sample by our observations (~ 194 arcmin 2 ); 
this corresponds to 0.026 spurious HB stars per arcmin 2 . On the basis of this, 11 field stars 



The large dispersion in the redder HB stars arises because RR Lyrae variables are included. 
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are expected to "contaminate" the HB population over the sampled cluster region (r < r t ). 



4. THE BSS RADIAL DISTRIBUTION 

The radial distribution of BSS in M5 has been studied following the same procedure 
previously adopted for other clusters (see references in Ferraro 2006; Beccari et al. 2006). 

First, we have compared the BSS cumulative radial distribution to that of HB stars. A 
Kolmogorov-Smirnov test gives a ~ 10~ 4 probability that they are extracted from the same 
population (see Figure [9]). BSS are more centrally concentrated than HB stars at ~ 4<r level. 

For a more quantitative analysis, the surveyed area has been divided into 8 concentric 
annuli, with radii listed in Table [2j The number of BSS (iVess) and HB stars (Nub), as 
well as the fraction of sampled luminosity (£ sam P) have been measured in the 8 annuli and 
the obtained values are listed in Table [2J Note that HB star counts listed in the table are 
already decontaminated from field stars, according to the procedure described in Section [3721 
(1, 2, and 8 HB stars in the three outer annuli have been estimated to be field stars). The 
listed values have been used to compute the specific frequency F^g S = A^bss/-^hb, and the 
double normalized ratio (see Ferraro et al. 1993): 



j-, V P°P/ pop/ /-i \ 

-"-pop — / rsamn /r sam P N > \ L ) 



with pop = BSS, HB. 



In the present study luminosities have been calculated from the surface density profile 
shown in Figure HI The surface density has been transformed into luminosity by means of a 
normalization factor obtained by assuming that the value obtained in the core (r < 27") is 
equal to the sum of the luminosities of all the stars with V < 20 lying in this region. The 
distance modulus quo ted in Section [2741 and a reddening E(B — V) = 0.03 have been adopted 



([Ferraro et al.lll999bl ). The fraction of area sampled by the observations in each annulus has 
been carefully computed, and the sampled luminosity in each annulus has been corrected for 
incomplete spatial coverage (in the case of annuli 3 and 8; see Figures [JJ and [2]) . 

The resulting radial trend of -Rhb is essentially constant with a value close to unity 
over the surveyed area (see Figure [10]). This is just what expected on the basis of the 
stellar evolution theory, which predicts that the fraction of stars in any post-main sequence 
evolutionary stage is strictly proportional to the fraction of the sampled luminosity (Renzini 
& Fusi Pecci 1988). Conversely, BSS follow a completely different radial distribution. As 
shown in Figure [10] the specific frequency Rbss is highly peaked at the cluster center (a 
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factor of ~ 3 higher than J?hb in the innermost bin), decreases to a minimu ntf at r ~ 10 r r , 
and rises again outward. The same behavior is clearly visible also in Figure [TT], where the 
population ratio A^bss/^Yhb is plotted as a function of r/r c . 

Note that the region between 800" and r t ~ 1290" (and thus also BSS-59, that lies at 
r ~ 995"5) has not been considered in the analysis, since our observations provide a poor 
sampling of this annulus: only 35% of its area, corresponding to ~ 0.4% of the total sampled 
light, is covered by the WFI pointing. However, for sake of completeness, we have plotted in 
Figure d2] the corresponding value of F£$g even for this annulus (empty circle in the upper 
panel): as can be seen, there is a hint for a flattening of the BSS radial distribution in the 
cluster outskirts. 



4.1. Dynamical simulations 



Following the same approach as iMapelli et al.l (120041 . 2006), we now exploit dynamical 



simulations to derive some clues about the BSS formation mechanisms from their observed 
radial distribution. We use the Monte-Carl o simulation c o de or iginally developed by Sig- 



urdsson & Phinney (1995) and upgraded in IMapelli et al.l (12004! . 2006). In any simulation 



run we follow the dynamical evolution of iV BSS within a background cluster, taking into 
account the effects of both dynamical friction and distant encounters. We identify as COL- 
BSS those objects having initial positions £ r c , and as MT-BSS stars initially lying at 
fi ^> t c (this because stellar collisions are most probable in the central high-density regions 
of the cluster, while primordial binaries most likely evolve in isolation in the periphery). 
Within these two radial ranges, all initial positions are randomly generated following the 
probability distribution appropriate for a King model. The BSS initial velocities are ran- 
domly extracted from the cluster velocity distribution illustrated in Sigurdsson & Phinney 
(1995), and an additional natal kick is assigned to the COL-BSS in order to account for the 
recoil induced by the encounters. Each BSS has characteristic mass M and lifetime ti ast . We 
follow their dynamical evolution in the cluster (fixed) gravitational potential for a time £, 
(i = 1, N), where each ti is a randomly chosen fraction of ti ast . At the end of the simulation 
we register the final positions of BSS, and we compare their radial distribution with the 
observed one. We repeat the procedure until a reasonable agreement between the simulated 
and the observed distributions is reached; then, we infer the percentage of collisional and 
mass-transfer BSS from the distribution of the adopted initial positions in the simulation. 

For a detailed discussion of the ranges of values appropriate for these quantities and 



7 Note that no BSS have been found between 3/5 and 5'. 
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their effects on the final results we refer to iMapelli et al.l (120061 ). Here we only list the 
assumptions made in the present study: 

- the background cluster is approximated with a multi-mass King model, determined as 
the best fit to the observed profil^]. The cluster central velocity dispersion is set to 
a = 6.5 km s -1 (Dubath et al. 1997), and, assuming 0.5M Q as the average mass of the 
cluster stars, the central stellar density is n c = 2 x 10 4 pc~ 3 (Pryor & Meylan 1993); 

- the COL-BSS are distributed with initial positions r» < r c and are given a natal kick 
velocity of 1 x a; 

- initial positions ranging between 5 r c and r t (with the tidal radius r t ~ 48 r c ) have been 
considered for MT-BSS in different runs; 



- BSS masses have been fixed to M = 1.2 M & (IFerraro et al.ll2006al ). and their charac- 
teristic lifetime to ti ast = 2 Gyr; 

- in each simulation run we have followed the evolution of N = 10, 000 BSS. 

The simulated radial distribution that best reproduces the observed one (with a reduced 
X 2 — 0.6) is shown in Figure [11] (solid line) and is obtained by assuming that ~ 80% of the 
BSS population was formed in the core through stellar collisions, while only ~ 20% is made 
of MT-BSS. A higher fraction ( £ 40%) of MT-BSS does not correctly reproduce the steep 
decrease of the distribution and seriously overpredict the number of BSS at r ~ 10 r c , where 
no BSS at all are found, but it nicely matches the observed upturning point at r ~ 13 r c (see 
the dashed line in Figure [Til. On the other hand, a population of only COL-BSS is unable to 
properly reproduce the external upturn of the distribution (see the dotted line in Figure [TT]) . 
and 100% of MT-BSS is also totally excluded. Assuming heavier BSS (up to M = 1.5 M Q ) or 
different lifetimes ti ast (between 1 and 4 Gyr) does not significantly change these conclusions, 
since both these parameters mainly affect the external part of the simulated BSS distribution. 
Thus, an appreciable effect can be seen only in the case of a relevant upturn, and negligible 
variations are found in the best-fit case and when assuming 100% COL-BSS. The effect starts 
to be relevant in the simulations with 40% or more MT-BSS, which are however inconsistent 
with the observations at intermediate radii (see above). 



By using the simulations and the dynamical friction timescale (from, e.g.. IMapelli et al. 



20061 ). we have also computed the radius of avoidance of M5. This is defined as the char- 



acteristic radial distance within which all MT-BSS are expected to have already sunk to 



8 By adopting the same mass groups as those of IMapelli et al.l (|2006l ). the resulting value of the King 
dimensionless central potential is Wo = 9.7 
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the cluster core, because of mass segregation processes. Assuming 12 Gyr for the age of 
M5 (Sandquist et al. 1996) and 1.2 M for the BSS mass, we find that r avoid ~ 10 r c . This 
nicely corresponds to the position of the minimum in the observed BSS radial distribution, 
in agreement with the findings of Mapelli et al. (2004, 2006). 

5. SUMMARY AND DISCUSSION 

In this paper we have used a combination of HST UV and optical images of the cluster 
center and wide-field ground-based observations covering the entire cluster extension to de- 
rive the main structural parameters and to study the BSS population of the galactic globular 
cluster M5. 

The accurate determination of the cluster center of gravity from the high-resolution 
data gives a(J2000) = 15 h 18 m 33!53, <5(J2000) = +2°4'57"06, with a la uncertainty of 0"5 
in both a and S. The cluster density profile, determined from direct star counts, is well fit 
by a King model with core radius r c = 27" and concentration c = 1.68, thus suggesting that 
M5 has not yet suffered the core collapse. 

The BSS population of M5 amounts to a total of 59 objects, with a quite asymmetric 
projected distribution (see Figure [Q and a high degree of segregation in the cluster center. 
With respect to the sampled luminosity and to HB stars, the BSS radial distribution is 
bimodal: highly peaked at r £ r c , decreasing to a minimum at r ~ 10 r c , and rising again 
outward (see Figures [TOl and ITT1) . 

The comparison with results of W06 has revealed that 43 (out of 59) bright BSS iden- 
tified by these authors at r ^ 450" are in common with our sample. Moreover, 4 additional 
stars classified as faint BSS in their study are in common with our BSS sample at r ^ 450". 
Considering that we find 56 BSS within the same radial distance from the center, this corre- 
sponds to 84% matching of our catalogue. The discrepancies are explained by different data 
reduction procedures, photometric analysis, and adopted selection criteria, other than the 
spurious identification of 8 BSS by W06, due a strongly saturated star in their sample. The 
central peak of the -Rbss distribution in our study is slightly higher (but compatible within 
the error bar) compared to that of W06, and we extend the analysis to larger distance from 
the center (out to r > 800"), thus unveiling the external upturn and the possible flattening 
of the BSS distribution in the cluster outskirts. 

Moreover, we have compared the BSS radial distribution of M5 with that observed 
in other GCs studied in a similar way. In Figure [12] we plot the specific frequency F^<? s 
as a function of (r/r c ) for M5, M3, 47 Tuc, and NGC 6752. Such a comparison shows 
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that the BSS radial distributions in these clusters are only qualitatively similar, with a 
high concentration at the center and an upturn outward. However, significant quantitative 
differences are apparent: (1) the F^ s peak value, (2) the steepness of the decreasing branch of 
the distribution, (3) the radial position of the minimum (marked by arrows in the figure), and 
(4) the extension of the "zone of avoidance," i.e., the intermediate region poorly populated 
by BSS. In particular M5 shows the smallest F^ s peak value: it turns out to be ~ 0.24, 
versus a typical value ^ 0.4 in all the other cases. It also shows the mildest decreasing slope: 
at r « 2 r c the specific frequency in M5 is about a half of the peak value, while it decreases 
by a factor of 4 in all the other clusters. Conversely, it is interesting to note that the value 
reached by F^<? s in the external regions is ~ 50-60% of the central peak in all the studied 
clusters. Another difference between M5 and the other systems concerns the ratio between 
the radius of avoidance and the tidal radius: r avoi d — 0.2 r t for M5, while r avoi( j ~ 0.13 r t for 
47 Tuc, M3, and NGC 6752 (see Tables 1 and 2 in Mapelli et al. 2006). 

The dynamical simulations discussed in Section H~TJ suggest that the majority of BSS in 
M5 are collisional, with a content of MT-BSS ranging between 20% and 40% of the overall 
population. This fraction seems to be smaller than that (40-50%) derived for M3, 47 Tuc 



and NGC 6752 by iMapelli et al.l (120061 ). in qualitative agreement with the smaller value of 
r avoid /r 4 estimated for M5, which indicates that the fraction of cluster currently depopulated 
of BSS is larger in this system than in the other cases. More in general, the results shown 
in Figure [TT] exclude a pure collisional BSS content for M5. 

Our study has also revealed the presence of a candidate BSS at ~ 24' from the center, 
i.e., beyond the cluster tidal radius (see Figures [2] and [6] and BSS-59 in Tabled]). If confirmed, 
this could represent a very interesting case of a BSS previously belonging to M5 and then 
evaporating from the cluster (a BSS kicked off from the core the because of dynamical 
interactions?). 
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Fig. 1— Map of the HST sample. The heavy solid line delimits the HST-WFPC2 FoV 
of our UV observations (Prop. 6607), the dashed line bounds the FoV of the optical HST- 
WFPC2 observations by Piotto et al. (2002), and the dotted line marks the edge of the 
complementary ACS data set. The derived center of gravity C grav is marked with a cross. 
BSS (heavy dots) and the concentric annuli used to study their radial distribution (cfr. Table 
[5]) are also shown. The inner and outer annuli correspond to r = r c = 27" and r = 115", 
respectively. 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Fig. 2. — Map of the WFI sample. All BSS detected in the WFI sample are marked as heavy 
dots, and the concentric annuli used to study their radial distribution are shown as solid lines, 
with the inner and outer annuli corresponding to r = 115" and r = 800", respectively (cfr. 
Table [2]). The circle corresponding to the tidal radius (r t ~ 21(5) is also shown as dashed- 
dotted line. The BSS lying beyond r t might represent a BSS previously belonging to M5 
and now evaporating from the cluster. 
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Fig. 3.— Optical CMDs of the WFPC2-HST and the WFI samples. The hatched regions 
indicate the magnitude limit (V < 20) adopted for selecting the stars used to construct the 
cluster surface density profile. 
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Fig. 4. — Observed surface density profile (dots and error bars) and best-fit King model (solid 
line). The radial profile is in units of number of stars per square arcseconds. The dotted line 
indicates the adopted level of the background, and the model characteristic parameters (core 
radius r c , concentration c, dimensionless central potential W^o) are marked in the figure. The 
lower panel shows the residuals between the observations and the fitted profile at each radial 
coordinate. 
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Fig. 5. — CMD of the ultraviolet HST sample. The adopted magnitude limit and selection 
box used for the definition of the BSS population are shown. The resulting fiducial BSS are 
marked with emp t y circ les. The open square corresponds to the variable BSS identified by 



Drissen fc Shara I (119981 ). The box adopted for the selection of HB stars is also shown. 
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Fig. 6.— CMD of the optical HST-WFPC2 and WFI samples. The adopted BSS and HB 
selection boxes are shown, and all the BSS identified in these samples are marked with the 
empty circles. The two BSS not included in the box in the left-hand panel lie well within 
the selection box in the UV plane and are therefore considered as fiducial BSS. The empty 
triangle in the right-hand panel corresponds to the BSS identified beyond the cluster tidal 
radius, at r ~ 24'. 
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Fig. 7. — CMD of the ACS complementary sample. The BSS selection box is shown, and 
the resulting fiducial BSS are marked with empty circles. 
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Fig. 8. — Left-hand panel: position of the 8 false BSS (marked with white circles) as derived 
from Table 2 of W06, overp lotted to the CFHT image (units are the same as in their Figure 
1). As can be seen, a heavily saturated star is responsible for the false identification. Right- 
hand panel: location of the 8 false BSS (empty circles) in the (B, B — I) plane, as derived 
from Table 2 of W06 (cfr. to their Fig. 2). 
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Fig. 9. — Cumulative radial distribution of BSS (solid line) and HB stars (dashed line) as 
a function of the projected distance from the cluster center for the combined HST+WFI 
sample. The two distributions differ at ~ 4<r level. 
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Fig. 10. — Radial distribution of the BSS and HB double normalized ratios, as defined in 
equation (CQ), plotted as a function of the radial coordinate expressed in units of the core 
radius. -Rhb (with the size of the rectangles corresponding to the error bars computed 
as described in Sabbi et al. 2004) is almost constant around unity over the entire cluster 
extension, as expected for any normal, non-segregated cluster population. Instead, the radial 
trend of -Rbss (dots with error bars) is completely different: highly peaked in the center (a 
factor of ~ 3 higher than -Rhb), decreasing at intermediate radii, and rising again outward. 
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Fig. 11. — Observed radial distribution of the specific frequency N BSS /Nn B (filled circles 
with error bars), as a function of r/r c . The simulated distribution that best reproduces the 
observed one is shown as a solid line and is obtained by assuming 80% of COL-BSS and 20% 
of MT-BSS. The simulated distributions obtained by assuming 40% of MT-BSS (dashed line) 
and 100% COL-BSS (dotted line) are also shown. 
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Fig. 12. — Radial distribution of the population ratio Nbss/N-hb for M5, M3, 47 Tuc, and 
NGC 6752, plotted as a function of the radial distance from the cluster center, normalized to 
the core radius r c (from Mapelli et al. 2006, r c ~ 30", 21", 28" for M3, 47 Tuc, and NGC 6752, 
respectively). The arrows indicate the position of the minimum of the distribution in each 
case. The outermost point shown for M5 (empty circle) corresponds to BSS-58, lying at 
r ~ 995". This star has not been considered in the quantitative study of the BSS radial 
distribution since only a negligible fraction of the annuls between 800" and r t is sampled by 
our observations. 
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Table I. The BSS population of M5 



Name 


RA 


DEC 


^255 


U 


B 


V 


I W06 




[degree] 


[degree] 













BSS-1 

BSS-2 

BSS-3 

BSS-4 

BSS-5 

BSS-6 

BSS-7 

BSS-8 

BSS-9 

BSS-10 

BSS-11 

BSS-12 

BSS-13 

BSS-14 

BSS-15 

BSS-16 

BSS-17 

BSS-18 

BSS-19 

BSS-20 

BSS-21 

BSS-22 

BSS-23 

BSS-24 

BSS-25 

BSS-26 

BSS-27 

BSS-28 

BSS-29 

BSS-30 

BSS-31 

BSS-32 

BSS-33 

BSS-34 

BSS-35 

BSS-36 

BSS-37 

BSS-38 



229.6354506 
229.6388102 
229.6383433 
229.6416234 
229.6416518 
229.6381953 
229.6403657 
229.6412279 
229.6376256 
229.6401139 
229.6396566 
229.6432834 
229.6384406 
229.6274500 
229.6204246 
229.6209379 
229.6264834 
229.6368731 
229.6367309 
229.6345837 
229.6382677 
229.6340227 
229.6332685 
229.6366685 
229.6393544 
229.6378381 
229.6349851 
229.6397645 
229.6370495 
229.6358816 
229.6361653 
229.6339822 
229.6281392 
229.6241278 
229.6332759 
229.6270877 
229.6244175 
229.6180419 



2.0841090 
2.0849660 
2.0842640 
2.0851791 
2.0836794 
2.0810119 
2.0824062 
2.0823768 
2.0793288 
2.0794858 
2.0784944 
2.0797197 
2.0776614 
2.0864896 
2.0879629 
2.0917858 
2.0960870 
2.0896002 
2.0917639 
2.0906438 
2.0934706 
2.0853879 
2.0875294 
2.0807168 
2.0762832 
2.0779999 
2.0807202 
2.0736403 
2.0770798 
2.0747883 
2.0720147 
2.0723032 
2.0756490 
2.0750261 
2.0603761 
2.0662947 
2.0693612 
2.0724090 



16.52 
17.95 
18.21 
17.59 
16.28 
17.36 
17.40 
17.91 
17.84 
17.57 
17.51 
18.12 
17.36 
18.07 
18.33 
17.80 
16.32 
16.56 
18.27 
17.88 
18.25 
17.67 
17.69 
18.23 
18.11 
17.86 
18.17 
18.19 
16.83 
18.25 
18.29 
16.73 
17.74 
18.21 
17.48 
17.33 
16.89 
17.37 



16.15 
17.38 
17.63 
17.22 
15.99 
16.99 
17.07 
17.47 
17.12 
16.98 
17.20 
17.64 
16.88 
17.63 
17.61 
17.28 
16.22 
16.30 
17.35 
16.81 
17.58 
17.32 
17.34 
17.78 
17.79 
17.52 
17.51 
17.60 
16.56 
17.81 
17.77 
16.10 
17.41 
17.50 
17.17 
17.18 
16.41 
17.23 



15.88 
17.40 
17.64 
17.05 
15.79 
16.81 
16.97 
17.41 
16.99 
16.87 
17.12 
17.78 
16.88 
17.64 
17.75 
17.26 
16.20 
16.11 
17.58 
16.96 
17.71 
17.22 
17.21 
17.67 
17.72 
17.43 
17.74 
17.69 
16.57 
17.79 
17.81 
16.16 
17.22 
17.65 
16.95 
17.06 
16.51 
17.12 



15.71 
17.04 
17.32 
16.90 
15.70 
16.65 
16.76 
17.15 
16.77 
16.62 
16.92 
17.54 
16.59 
17.33 
17.36 
16.98 
16.13 
16.01 
17.07 
16.43 
17.35 
17.03 
17.08 
17.37 
17.50 
17.27 
17.30 
17.28 
17.75 
17.51 
17.47 
15.95 
17.09 
17.27 
16.86 
16.95 
15.71 
17.00 



CR2 

CR4 

CR3 

CR5 

CR1 

CR21 

CR12 

CR13 

CR23 

CR22 

CR24 

CR25 

CR8 

CR11 

CR18 

CR20 

CR14 

CR17 

CR16 

CR19 

CR7 

CR10 



CR70 
CR33 
CR26 
CR31 
CR36 
CR35 
CR29 
CR79 
CR48 
CR47 
CR46 
CR34 



-31 - 



BSS- 


-39 


229. 


.6311963 


2 


BSS- 


-40 


229. 


.6297499 


2 


BSS- 


-41 


229. 


.6443367 


2 


BSS- 


-42 


229. 


,6448646 


2 


BSS- 


-43 


229. 


,6460645 


2 


BSS- 


A A 

-44 


229. 


,6481631 


2 


BSS- 


A V 

-45 


229. 


,6433942 


2 


BSS- 


-46 


229. 


.6439884 


2 


BSS- 


-47 


229 


.6180420 


2 


BSS- 


-48 


229. 


.6092873 


2 


BSS- 


-49 


229 


.6723094 


2 


BSS- 


-50 


229 


,6006551 


2 


BSS- 


-51 


229 


,6669956 


1 


BSS- 


-52 


229 


.5949935 


2 


BSS- 


-53 


229 


,6706625 


2 


BSS- 


-54 


229 


,6667908 


2 


BSS- 


-55 


229 


.7370667 


2 


BSS- 


-56 


229 


.5476990 


2 


BSS- 


-57 


229 


.6711255 


1 


BSS- 


-58 


229 


.4381714 


2 


BSS- 


-59 


229 


.7408412 


2 


BSS- 


-60 


229. 


.3218200 


2 



Table 1 — Continued 

0857800 18.31 17.33 17 

0664961 18.16 17.58 

0872809 17 

0738335 - 16. 

0748695 - - 16 

0718829 - - 16. 

0760163 - - 17 

0775670 17. 

0598328 

1680914 16 

0882827 16. 

0814678 - - 17 

9781808 - - 17. 

0469325 - - 17. 

0695464 - - 17. 

1149550 - 17. 

0323392 17 

0112610 - - 16 

9415566 16 

0302088 - - 17 

3399166 17 

3271022 - - 16 



40 


16.76 








17.27 




CR76 


50 


17.23 


— 


CR9 


53 


16.06 


15.95 


CR32 


64 


16.44 


16.66 


CR30 


72 


16.61 


16.87 


CR37 


03 


16.79 


16.91 


CR28 


44 


16.99 


16.81 


— 




17.18 


17.12 




85 


16.68 




0R2 


94 


16.64 




0R9 


00 


16.74 




OR10 


20 


16.74 




0R1 


69 


17.46 




0R4 


82 


17.50 






82 


17.72 






80 


17.42 




OR23 


88 


16.60 




0R5 


98 


16.64 






75 


17.33 






49 


17.08 






34 


16.09 







Note. - - The first 41 BSS have been identified in the WFPC2 sample; BSS-42- 
46 are from the complementary ACS observations; BSS-47-59 are from the WFI 
data-set. BSS-59 lies beyond the cluster tidal radius, at ~ 24' from the center. The 
last column list the corresponding BSS in W06 sample, with "CR" indicating their 
"Core BSS" and "OR" their "Outer Region BSS". 
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Table 2. Number counts of BSS and HB 
stars 



„.// 


r e " 


^BSS 


N H b 


L samp /L: 





27 


22 


94 


0.14 


27 


50 


15 


94 


0.16 


50 


115 


10 


135 


0.26 


115 


150 


3 


46 


0.09 


150 


210 


2 


52 


0.10 


210 


300 





45 f 


0.10 


300 


450 


4 


42* 


0.09 


450 


800 


2 


38 f 


0.06 



Note. - t The A^hb values listed here 
are those corrected for field contamination 
(i.e., 1, 2 and 8 stars have been subtracted 
to the observed number counts in these 
three external annuli, respectively). 



